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ABSTRACT 


The purpose of this document is to present the results of 
comparisons of the solar flux models. (The wavelength X = 
10.7 cm radio flux is the best indicator of the strength of 
the ionizing radiations such as solar ultraviolet and x-ray- 
emissions that directly affect the atmospheric density 
thereby changing the orbit lifetime of satellites. Thus, 
accurate forecasting of solar flux F^q 7 is crucial for 
orbit determination of spacecrafts.) The measured solar 
flux recorded by National Oceanic and Atmospheric Admini- 
stration (NOAA) is compared against the forecasts made by 
Schatten, Marshall Space Flight Center (MSFC) , and NOAA it- 
self. This document also discusses the possibility of a 
combined linear, unbiased minimum-variance estimation that 
properly combines all three models into one that minimizes 
the variance. All the physics inherent in each model are 
combined. This is considered to be the dead-end statistical 
approach to solar flux forecasting before any nonlinear cha- 
otic approach. 

The research for this document was completed in December 1990. 
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SECTION 1 - INTRODUCTION 


This document is the first part of a sequence of preliminary- 
studies of solar flux observed at the wavelength X = 

10.7 cm range. The sequence starts with comparisons of dif- 
ferent solar flux models and gradually leads to a critical 
stochastic approach, which further produces a geometric 
approach to the prediction of chaotic solar flux time series. 

The analysis in this first sequence is based on the avail- 
able forecasts by Schatten (at Goddard Space Flight Center 
(GSFC) ) (Reference 1), Marshall Space Flight Center (MSFC) 
(Reference 2), and National Oceanic and Atmospheric Admini- 
stration (NOAA) (Reference 3). The comparisons are made 
against actual observed values that are collected by NOAA 
(Reference 3 ) . 

The observable radio spectrum extends from 1 centimeter (cm) 
to 10 meters (m) . Like the optical spectrum, the radio 
spectrum is limited on its short wavelength end by absorp- 
tion in the Earth's atmosphere (by molecules of oxygen and 
water vapor) . On the long wavelength end, the lower atmos- 
phere is always transparent, even on cloudy days. But a 
high layer, called the ionosphere, begins to interfere at 
around X = 10 m (References 4 and 5). 

The radio waves are radiated by fast-moving electrons in the 
highly ionized gases of the outer solar atmosphere. Ionized 
gases, which are fully transparent to visible light, how- 
ever, may be opaque to radio waves at certain wavelengths. 

The opacity depends on the density of ionized gas. In the 
solar chromosphere, where density is high, the gases are 
completely opaque to meter wavelengths; only the centimeter 
waves can escape the Sun to reach the Earth. The Sun that 
is observed is only the visible Sun; it appears larger in 
the radio region (that is, the appearance dimension is pro- 
portional to the wavelength) . 
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There is a strong correlation between sunspots and the solar 
flux F 1q 7 because probably most of the enhanced radiation 
comes from limited areas of the Sun where there are active 
sunspots. The activity depends on wavelength of radiated 
solar flux. For waves shorter than 3 cm, the intensity is 
steady. From 3 to 60 cm, often called decimeter range, the 
intensity shows occasional short-lived increases. These 
tend to last for a few minutes. The decimeter intensity 
also shows a slowly varying component that tends to exhibit 
a 27-day period associated with solar rotation (Reference 6) 
and rises from the vicinity of active sunspot regions. 

Large sporadic outbursts, lasting for minutes, occur often 
in association with the bigger solar flares. A millionfold 
increase in intensity within a few seconds has been observed. 
(This will be studied as a part of the sequence of the solar 
flux analysis by identifying the abrupt changes as one of 
seven Thom's "elementary catastrophes.") (See Section 5 for 
recommendations . ) 

It is necessary to study solar flux and accurately forecast 
it to perform accurate orbit determination for a spacecraft. 
The orbit lifetime is a function of atmospheric drag force; 
this force is a function of atmospheric density, which it- 
self is a function of solar flux. 

Section 2 is devoted mostly to graphical analysis of solar 
flux data. Section 3 describes the statistical techniques 
to compare different forecasting models by confidence inter- 
val methods. Section 4 introduces a linear, unbiased 
minimum-variance estimation and combines three important 
models into one. Section 5 is devoted to conclusions and 
recommendations for future investigations. 
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SECTION 2 


SOLAR FLUX F 


PREDICTION 


10.7 


2.1 MATHEMATICAL PERSPECTIVES 

satellite orbit lifetime prediction, one has to evaluate 
the drag force that continually results in satellite orbit 
decay. By applying the fundamental principles of fluid 
mechanics, the drag force is written in the following form: 

|F| = ^ p V 2 C a A (2-1) 

where p is the atmospheric density, which is a complicated 
function of solar flux in different density models of the 
atmosphere. The velocity of the spacecraft is indicated by 
v ' and th e other variables are properties of the space- 
craft. These properties are drag coefficient C d and scat- 
tering cross section A. 

It is very clear from the above equation that, given the 
drag coefficient and the scattering cross section A, 
one can easily calculate the drag force fF| if the den- 
sity of the atmosphere is known. Since the atmospheric den- 
sity is sensitive to solar activity, most of the density 
models are complicated functions of solar flux. The motiva- 
tion for studying solar flux prediction models (other than 
solar astronomy) is the accurate satellite orbit lifetime 
prediction. 

2 - 2 GRAPHICAL ANALYSIS OF DATA 

The first part of this study compares Schatten solar flux 
forecasts with National Oceanic and Atmospheric Administra- 
tion (NOAA) actual solar flux values. Schatten predictions 
are modified inconsistently, but at least once every 

3 months. NOAA also forecasts short-term predictions that 
are modified consistently every week. Schatten’ s 
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latest values were distributed on September 1, 1990, which 
includes long-term predictions starting from September 1990 
to September 2008. His next latest predictions were dis- 
tributed on May 25, 1990, which includes predictions start- 
ing May 1990 to August 1990 and on to April 2012. Therefore 
his May 25 version, which includes May, June, July, and 
August 1990 predictions, is by far the best he could do. 

Thus this analysis was done on his best predictions (May, 
June, July, and August 1990). 

As seen from the graphs in Figures 2-1 to 2-10, the 30-day 
mean of the actual (NOAA) values are always less than 
+2 sigma value of Schatten and in most cases are even 
smaller than Schatten (mean) . They are closer to -2 sigma 
value or even smaller than that. This may mean that the 
mission analysis is being too conservative by using +2 sigma 
value consistently. 

2.3 DESCRIPTION OF GRAPHS 

Figures 2-1 to 2-10 present the actual solar flux values and 
different forecasts for the months of May, June, July, and 
August 1990. The forecasts are the best updates for those 
months. Figures 2-11 to 2-14 present the actual solar flux 
values and their different averages and different forecasts 
for a period of 2 years (October 15, 1988 to September 17, 
1990). The actual data are daily values, and the long-term 
prediction models (Schatten and MSFC) are monthly values. 
Figures 2-15 and 2-16 are the confidence intervals for 
Schatten and MSFC forecasts, the nominal and the +2 a, re- 
spectively. Figures 2-17 to 2-19 are the actual solar flux 
values and their averages for three different timespans. 

The statistical analysis performed to get the confidence 
intervals is discussed in Section 3. It should be noted 
that all the units for the solar flux values are in units of 
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actual(NOAA) 


Data from "Science Data" * 



0 10 20 30 40 


MAY 1990 


<3 actual (NOAA) 

• Schatten mean 

Schatten -2 a 

Schatten +2 a 


* Science data is a data file with its data printed by 
running a cricket graph on the Macintosh. 


Figure 2-1. Plot of Solar Flux Values and Schatten 
Nominal, +2 Sigma, and -2 Sigma Pre- 
dictions (May 1990) 
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Data from "solar flux comparison 
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Data from "solar flux comparison 
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Figure 2-3. Plot of Solar Flux Values, Schatten Nominal and -2 Sigma 

Predictions, and NOAA 27-Day Nominal Predictions (May 1990) 



act ual-mean/slg 


Data from "solar flux comparison" 



Figure 2-4. Plot of Solar Flux Values Minus the Mean 
Value Normalized by Standard Deviation 
(May 1990) 
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Data from "Science Data 
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Figure 2-6. Plot of Solar Flux Values and the Schatten Nominal, -2 Sigma 
Predictions, and the NOAA 27-Day Nominal Predictions 
(June 1990) 
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Sigma Predictions (July 1990) 



Data from "Science Data 
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Figure 2 8. Plot of Solar Flux Values, the Schatten Nominal, -2 Sigma 
Predictions, and the NOAA 27-Day Nominal Predictions 
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+2 Sigma Predictions (August 1990) 
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Figure 2 10. Plot of Solar Flux Values, the Schatten Nominal, -2 Sigma, 
and NOAA 27-Day Nominal Predictions (August 1990) 




Figure 2-11. Plot of Actual NOAA Solar Flux Values for October 15 
to September 17, 1990 



solar flux 



Figure 2-12. Plot of Actual Solar Flux Values and Mean Value 







3001 


i 


CM 

0 $ 


o 

CT> 



ssrifBA xnyjeios 


6176 


2-15 


Figure 2 13. Plot of Solar Flux Values and the Nominal Best Schatten and 
MSFC Predictions 




Actual (NOAA) Best Schatten (NOM) Best MSFC (NOM) 

x = 200.04 x = 181.75 " a " x = 197.27 


Best MSFC (97.7%) 
x = 212.56 


Best Schatten (+2) 
x = 218.85 


Figure 2-14. Plot of Solar Flux Values, the Schatten 
Nominal and +2 Sigma Predictions, and 
the MSFC Nominal and 97.7 Percent 
Predictions 
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Confidence Interval 


Schatten Nominal Confidence lnterval(95%) 
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LCI(Sch-NOM) = 170.03 _ 
UCI(Sch-NOM) = 193.47 X = 181 - 75 
Actual Mean = 200.04 


MSFC Nominal Confidence Interval (95 
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Figure 2-15. Plot of the Nominal Schatten and MSFC 
95-Percent Predictions Confidence 
Interval and the Average of the Actual 
Solar Flux Values 
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Confidence Interval 


Schatten(+2Sigm) Confidence lnterval(95%) 
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MSFC 97.7 Confidence Interval (95%) 
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Figure 2-16. Plot of the Schatten +2 Sigma and MSFC 
97.7-Percent Predictions Confidence 
Interval and the Average of the Actual 
Solar Flux Values 
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Figure 2 18. Plot of the Actual Solar Flux Values for a Different 
Timespan (MJD 13400 to 15200) 
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Figure 2 19. Plot of the Actual Solar Flux Values and Its Average for 
Different Timespan (MJD 15000 to 16800) 



2 -22 

watt/m /Hz x 10 for the range X = 10.7 cm wavelength, and 
the horizontal axis is time in modified Julian date. 

In order to compare forecast models of solar flux, one can 
compare the forecasts of solar flux (F 1Q 7 ) made for the 
timespan that the actual solar flux values are available. 
Every forecast will result in an interval (with a certain 
percentage of confidence) . The actual population mean of 
data will fall within that confidence interval. This inter- 
val can be calculated for each forecast model. To compare 
forecast models, the question is whether the confidence in- 
terval encloses the population mean of the actual solar flux 
values or not. If it does, 95.5 percent of the time the 
predicted value is within the confidence interval; thus, it 
is a good forecast. The mathematical analysis of this pro- 
cedure is presented in the next sections. 
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SECTION 3 - STATISTICAL ANALYSIS 


3.1 CHECK OF A HYPOTHESIS 

Sample values are often used as estimators for parameters of 
random variables. However, these procedures result only in 
point estimates for a parameter of interest; no indication 
is provided about how closely a sample value estimates the 
parameter. A more meaningful procedure for estimating 
parameters of random variables involves the estimation of an 
interval, as opposed to a single point value, which will 
include the parameter being estimated with a known degree of 
uncertainty. For example, consider the case where the sam- 
ple mean x computed from N independent observations of a 
random variable x is being used as an estimator for the mean 
value p x . It is usually more desirable to estimate p 
in terms of some interval d, such as x + d, where there is a 
specified uncertainty that p x falls within that interval. 
Such intervals can be established if the sampling distribu- 
tion of the estimator in question is known (Reference 7). 

It can be shown that probability statements can be made con- 
cerning the value of a sample mean x as follows. 


prob 


U - u x ) v/n 

Z l- (a/2) < o x - Z a/2 


1 - a (3-1) 


prob 


x > 



a 


(3-2) 


where x = sample mean, p x = population mean, N = number 
of observations, ct x = sample standard deviation, d = 
uncertainty length, a = probability measure, and z 

a 

desired percentage of confidence. 
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As the sample size N becomes large, the sampling distribu- 
tion of the sample mean x approaches a normal distribution 
regardless of the distribution of the original variable x. 

For a sample, the probability statement would be either 1 or 
0 , i . e . , 


(x - u ) v/n 

prob z. , ... < < z 

1 1- (a/2) a - a/2 


0 

u 


( 3 - 3 ) 


As the value a becomes small (as the interval between 

z l-(a/2) and z a/2 becomes wide), the probability is 
more likely to be unity rather than zero. In slightly dif- 
ferent terms, if many different samples were repeatedly col- 
lected and a value x were computed for each sample, one 
would expect 


prob 


0 

1 


( 3 - 4 ) 


to fall within the noted interval for about 1-a of the 
samples. In this context a statement can be made about an 
interval within which one would expect to find the quantity 


(x - p x ) yi 


( 3 - 5 ) 


with a small degree of uncertainty. Such a statement is 
called a confidence statement. The interval associated with 
the confidence statement is called a confidence interval. 
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For the case of the mean value estimate, a confidence inter- 
val can be established for the mean value p x based upon 
the sample value x by rearranging terms in the previous 
equation as follows: 


a x z a/2 . . — 

x - < u < x + 

— K x 


a z 
x a/2 




n/n 


(3-6) 


3.2 CHI-SQUARE DISTRIBUTION 

2 9 

For variance a based on the sample variance s for 
a sample size of N = 31, one would use chi-square distribu- 
tion x 2 


2 9 

n s 2 _ ns 

2 < a x < 2 

n;a/2 x n;l-a/2 


n = N 


1 


(3-7) 


30s 2 2 30s 2 

2 °x 2 

x 30;a/2 x n;l-a/2 


n = 31 - 1 = 30 


(3-8) 


From a standard statistical table called, “Percentage Points 
of Chi-Square Distribution" (Reference 8) the value of 
x n*a/2 Given a can be found. 

For the value of a = 0.10, 1 - a/2 = 0.95, and a/2 = 0.05, 

X 30;a/2 = 43 ' 77 and X 30;l-a/2 = 18 - 49 - 
So the interval reduces to 


[0.6854s 2 < cr 2 < 1.622s 2 ] (3-9) 

Calculate the sample mean x and the sample variance s 2 to 
find both intervals. 
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SECT ION 4 - LINEAR. UNBIASED MINIMUM- VARIANCE 
ESTIMATION ( LUMVE ^ 


Schatten prsdictions are now adopted by GSFC bscauss fchsy 
apparently did a good job at some periods of time. But the 
conclusion from the data for the past 2 years is that MSFC 
predictions were closer to the actual solar flux values. 

This conclusion shows that it is not possible to compare the 
accuracy of two forecasting models (which use stochastic 
methods) when they try to model a time series that is inher- 
ently chaotic (existence of a structure in data). Therefore 
the best method is to combine all the models into one. This 
method is investigated in this section under linear, un- 
biased minimum-variance estimation (LUMVE) . In this method, 
the three solar forecasting models — NOAA , MSFC, and Schatten 
predictions are combined into one that minimizes the vari- 
ance . 

4 . 1 MATHEMATICAL FORMULATION 

The LUMVE ensures that the variance of the combined solar 
flux predictions is the smallest that can be achieved for 
any linear, unbiased combination of the individual predic- 
tions. This method was used because it was demonstrated 
that MSFC predictions were closer to the actual solar flux 
than Schatten' s for the past 2 years. 

Solar flux is an inherently unpredictable phenomenon, and 
stochastic methods used by Schatten, MSFC, and NOAA cannot 
produce good predictions. Therefore, the dead-end approach 
before an analytic (not stochastic) nonlinear, chaotic ap- 
proach is the linear, unbiased minimum-variance approach 
(References 9 and 10). 

Let 


o n = standard deviation of NOAA prediction 
a m = standard deviation of MSFC prediction 
a s = standard deviation of Schatten prediction 
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Consider a linear combination of the three flux predictions. 
Using normalized coefficients 


f 


b f + b f + b f 
n n in m s s 

b + b + b 
n m s 


( 4 - 1 ) 


where b n , b m , and b s are coefficients for NOAA, 

and Schatten forecast of solar flux f , f . and 

n m 

respectively. 


MSFC , 


Define 


a 

n 


n 

b + b + b 
n m s 


( 4 - 2 ) 


a m = 


m 


b + b + b 
n m s 


( 4 - 3 ) 


a_ = 


b + b + b_ 
n m s 


( 4 - 4 ) 


thus 


/X 

f = 


af + 
n n 


a m f_ + 
m m 


a s f s 


( 4 - 5 ) 


Now the problem is to select values of a , a , and a 

^ n m s 

that will yield the best f (the closest value to the mean of 

the actual data) . 

/\ 

The variance of f is 


v(f) 


2 2 
a a 
n n 


,2 2 
a a + 
m m 


2 2 
a a 
s s 


+ covariance terms 
0 


for independent predic- 
tions 


( 4 - 6 ) 
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By imposing E(f) = f (where f is the desired actual mean), 
and forcing the parameters 


E ( f ) = E(f n ) = E(f m ) = E(f s ) = f (4-7) 


which is a renormalization technique. Then, 


E< £ > = a n E(£ n ) + a m E(£ m ) + a s E(f s ) (4-8) 


and 


f = a f + a f + a f 
n m s 


(4-9) 


1 — 3_ + 3 + a 

n m s 


(4-10) 


Equation (4-6) then becomes 


V(f) = a 


2 2 2 2 /, \ 
_ o + a or + [1 - a - a | 

n n m m I n ml 


a ,2 2 

- a | ct 
m | s 


(4-11) 


or 


V(f) 


2 2 2 \ 2 2 2 \ 

= a a + a + a m o m + a - 
n l n sj m l m s / 


2 2 

2a ct - 2a a 
ns ms 


2 2 

+ 2a a m a + u 
n m s s 


(4-12) 


If the variance V(f) is minimized (the ideal case), that is 


9V(f) 

3a 

n 


= 0 


(4-13) 


then * o*j a a m a* = and 
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9V(f) 

9a 

m 


= 0 


(4-14) 


then a n (o_) 2 + a m ( a m + ) = c ^ ; 

ti\ s ) m \ m s / s 


therefore 


a n = 


2 2 
a a 
m s 


2 2 2 2 2 2 

aa+aa+aa 
n m ns ms 


(4-15 ) 


a m = 


2 2 

cr 

n s 


2 2 2 2 2 2 

a^ + a 

n m ns ms 


(4-16) 


a 


2 2 
a a 
n m 


2 2 2 2 2 2 

a « G m + + 
n m ns ms 


(4-17) 


By dividing by a 2 a 2 a 2 
2 * 2 n m s 


1/a' 


a = 
n 


n 


0) 


2 2 2 
l/a„ + 1/at + l/a z 
n m s 


n 


a) +a) + aj 

n m s 


(4-18) 


1/a' 


a m " 


m 


Cl) 


I/O 2 + I/O 2 + I/O 2 

n m s 


m 


CO + (0 + CO 

n m s 


(4-19) 


and 


1 /a! 




a = 


1/aJ + l/a* + l/o 2 
n m s 


+ go 

n m s 


(4-20) 
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where 



f 


W f + 0) f + (0 f 
n n m m s s 

CO + 01 + 0) 

n m s 


(4-21) 


This approach requires only the ratios of the coefficients 
( a n' a m' a s ^ and not the actual parameters. See Fig- 
ure 4-1 for calculations of these parameters. The variances 
are calculated on the PC (IBM AT compatible) using the 
Quattro Pro program. 

Can the same coefficients be used in the future or do these 
coefficients vary with time? If they vary in time are the 
variations predictable or not? Figures 4-1 through 4-10 
show that these coefficients evolve in time in a predictable 
fashion. The reason they are predictable is that all the 
variations are already in the flux values and their adjust 
ment coefficients a^, a^, a^ do not vary violently and are 
predictable (Figures 4-1 through 4-10). 

4.2 GRAPHICAL ANALYSIS OF DATA 

The calculations of the required parameters in LUMVE pre- 
sented in the previous section is performed here. The inter- 
mediate coefficients — a n> a m , and a g — are calculated from 
Equations (4-15) through (4-21) and presented in Figure 4-1. 
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Figure 4 1. Calculation of the Parameters Used in Linear, Unbiased Minimum-Variance 
| Estimation (Running Procedure: Running Average of Data) 




Predicted Flux values 


300 


Predicted Flux 


ft ■ 


1 



Months from Jan 88 


NOAA-m 
MSFC-m 
Schatten-m 
Our Flux 


Time evolution of flux coefficients 



An 
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Figure 4-2, 


Predicted Flux by NOAA, MSFC, Schatten 
(Nominals) and the Result of Linear, 
Unbiased Minimum-Variance Estimation 
(Time Evolution of the Coefficients 
a n> a nt' an< 3 a s — Running Procedure) 
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Coefficients of Fiux 


Flux values predicted 
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Every six months from Jan 88 


Evolution of Flux coefficients 



Every six months from Jan 88 


Figure 4-4. 


Plot of the Predicted Flux by NOAA, MS 
Schatten (Nominal) and Result of Linea 
Unbiased Minimum-Variance Estimation 
(Time Evolution of the Coefficients 
a n* a m» an d a s — Partitioned 
Procedure) 
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Figure 4- 



Plot of the Actual Solar Flux Values, 
Schatten, and MSFC Predictions With 
30-Day Running Average and 30-Day 
Partitioned Average 
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Solar Flux Actual & Prediclion 
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Figure 4-9 


Plot of the Actual Solar Flux Values and the 
Schatten and MSFC Predictions With 30-Day 
Running Average and the Estimated Combined 
Flux With Partitioned Procedure 
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Finally the combined flux is calculated from Equation (4-21) . 
Note that these equations require variances of NOAA, MSFC, 
and Schatten data that had been calculated using the Quattro 
Pro program. 

The time evolutions of a n , a ffl , and a s are presented in 
Figures 4-2 through 4-4, which indicate that their varia- 
tions are not violent and are predictable in a sense. Fig- 
ures 4-5 through 4-8 present the actual solar flux values 
for a period of 2 years with MSFC and Schatten predictions 
and different kinds of averages, so that one can clearly see 
that MSFC predictions are closer to the actual values than 
Schatten predictions. Figure 4-9 presents the combined flux 
values, which are better than all the other models, because 
it apparently has all the inherent physics of the individual 

models built into it by the proper coefficients (a , a 
, . n' m' 

and a s ). Figure 4-10 presents the normalized deviations 

of the predictions by individual models from the 30-day 

average of the actual solar flux values. It is clear from 

this graph that the combined solar flux prediction model 

varies within 20 percent of the actual values, whereas the 

other individual prediction models show variations much 

larger than 20 percent from the actual values. 
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SECTION 5 - CONCLUSIONS AND RECOMMENDATIONS 


The statistical analysis using confidence intervals reveals 
that one cannot draw a general conclusion about which solar 
flux forecast model is better than the others among 
Schatten, NOAA, and MSFC. This is due to the fact that 
these models assume stochasticity (structurally random) in 
solar flux time series, which is chaotic (existence of an 
underlying structure in data). 

Before employing the nonlinear, chaotic approach that will 
follow as the second sequence of the analysis (to be pub- 
lished in a different document), a combined LUMVE has been 
developed that properly combines all three models into one 
that minimizes the variance. All the physics inherent in 
each model are combined. 

In the second part of these studies, solar flux as a chaotic 
time series will be studied and a particular route through 
which the dynamical time series becomes chaotic will be 
identified. The third part of the sequence is a critical 
stochastic approach to solar flux. In this part solar flux 
is studied through model identification, estimation, 
fitting, diagnostic checking, and mathematical forecasting 
for solar flux chaotic data. 

The following future investigations are recommended: 

• A Box-Jenkins type approach to solar flux time 
series model identification, estimation, fitting, diagnostic 
checking, and forecasting of solar flux. This is a method 
to classify the solar flux time series as one of the 
presently known models (moving average (MA) model, auto 
regressive (AR) model, and mixed auto regressive moving 
average (ARMA) model). Once the classification is made, 
proper forecastings seem possible. 
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• Chaotic approach to solar flux prediction. Since 
solar flux is shown to be chaotic, this approach would allow 
for the construction of an iterative manifold that can re- 
construct the solar flux time series. To do this a sequence 
of studies should be performed such as 

Finding Lyapunov spectrum from solar flux time 
series (Reference 11) 

Forming attractors from solar flux time series 
and nonlinear signal processing using Neural 
Net 

Extracting self similarity character and 
fractal structures from solar flux time series 
and modeling abrupt changes with Thom's ele- 
mentary catastrophes (Reference 12) 
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